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l/^ I We present self-similar models of resistive viscous Keplerian disks driv- 

ing non-relativistic magnetohydrodynamics (MHD) jets becoming super fast- 
J> ' magnetosonic. We show that in order to obtain such solutions, the thermal pres- 

ly-^ ' sure must be a sizeable fraction of the poloidal magnetic pressure at the Alfven 
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surface. These steady solutions which undergo a recoUimation shock causally 
disconnected from the driving engine, account for structures with a high temper- 
ature plasma in the sub-Alfvenic region. We suggest that only unsteady outflows 
with typical time-scales of several disk dynamical time-scales can be produced if 
Ph| the suitable pressure conditions are not fulfilled. 
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, 1. Introduction 

Self-coUimated jets are now commonly observed originating from young stellar objects 
(YSOs), active galactic nuclei and galactic binaries (Livio 1997). All these flows share 
common properties like being always correlated with the accretion phenomenon (Hartigan 
et al. 1995; Falcke & Biermann 1996; Serjeant et al. 1998; Gallo et al. 2003). It has long been 
identified that jet self-confinement requires the presence of a large scale magnetic field in order 
to focus the outflowing plasma (Chan & Henriksen 1980). The "universal" paradigm of jet 
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formation relies on the occurrence of bipolar magnetic fields threading the accretion disk. As 
a consequence, the theory of accretion disks had to be revisited in order to take into account 
the mass, angular momentum and energy extractions achieved by the jet. One notorious 
modification to the standard picture is the necessary radial stratification of the disk accretion 
rate, namely Ma oc being a measure of the disk ejection efficiency). For instance ^ = 
describes a standard disk with no outflow while < ^ < 1 stands for an ejecting Keplerian 
accretion disk (Ferreira 1997). If one wishes to obtain the exact ejection efficiency, one has 
to solve without any approximation the full magnetohydrodynamics (MHD) 2D structure of 
the disk. 

Anomalous magnetic diffusivity must be present whithin the disk to allow accreting 
(and rotating) mass to cross the magnetic field lines whereas ejected mass becomes frozen in 
to the field. Only self-similar solutions taking into account the underlying resistive accretion 
disk hitherto provided this description (Ferreira & Pelletier 1995; Ferreira 1997; Casse & 
Ferreira 2000a,b). Within these solutions, once in ideal MHD regime, mass is magnetically 
accelerated along each field line and must successively cross three MHD critical points, 
namely the slow and fast magnetosonic ones (SM,FM) and the Alfvenic point (A). So far, 
none of the self-similar solutions was able to obtain both disk and jet fiows, the latter crossing 
the three critical surfaces. 

Using the same framework, Vlahakis et al. (2000) solved the ideal MHD jet equations and 
provided new solutions crossing these three critical points. However since these solutions were 
not connected to the underlying disk, the issue of supcr-FM jet production from accretion 
disk remained. In this letter, we intend to present the necessary conditions to get self-similar 
super-FM jets (Sect. 2), discuss the properties of typical solutions (Sect. 3) and conclude 
with some astrophysical implications that may be put to the test of observations. 

2. Role of a sub-Alfvenic heating 

Stationary jets are described by a set of axisymmetric ideal MHD equations. Thus the 
poloidal magnetic field writes Bp = (Va x e^)/r, where a(r, z) = Cst describes a surface of 
constant magnetic fiux. Disk winds are produced whenever a large scale magnetic field, close 
to equipartition with the disk thermal pressure (Ferreira & Pelletier 1995), is present over a 
range in anchoring radii Tq. The corresponding jet is made of magnetic surfaces nested one 
around each other with several integrals of motion. In the non-relativistic case, one gets (up 
poloidal velocity, Q angular velocity and p density): (1) the mass to magnetic fiux ratio r){a) 
with Up = ?7(a)Bp/ /JoP', (2) the angular velocity of a magnetic surface Q^,{a) = Q — riB^/poP^ 
and (3) the specific total angular momentum L{a) = f2*r^ = flr^ — rB^/rj transported 



-3- 



away. Here, ta is the Alfven radius where mass reaches the Alfven poloidal velocity. In 
this letter, we are interested in jets that may be heated by their surroundings so that an 
adiabatic description is inadequate. Instead, we will assume the presence of a heat flux 
q = VH — VP/ p, where H is the usual enthalpy for a perfect gas. Including this additional 
effect, one gets the generalized Bernoulli invariant E{a) +J-'{s,a) = ^ + H + — f^^*B(f,/rj, 
where a) = q • e\\ds' is the heating term that depends on a curvilinear coordinate 
s along a given magnetic surface (s+ is roughly the SM point and Bp = -Bpe||). The total 
specific energy provided at the disk surface is E{a) ~ Qlrg{X — 3/2) for a thin disk, where 
is the Keplerian rotation at the anchoring radius To and A = L/Qor^ is the magnetic lever 
arm. The shape of the magnetic surface is given by the Grad-Shafranov (GS) equation 

(l-m2)J^ = Jx + J^ + Jf3, (1) 

f d£ ,^ 2*^^* dVt^r\ 
where Jx = pr \ — — I- (1 — gji^hr — h g\l 



^ da da da 

Bl-m^BldlnpA ^ 2^0 

= ^— — ; Vmp 

2/io da 

J0 = A(v^_q).e^ 

-Dp 

Here, = Va/|Va|, m = Up/VAp is the Alfvenic Mach number and g = 1 — Q/Q^. GS 
equation provides a(r, z) for a given set of invariants. Unfortunately, it is a PDE of mixed 
type: it is hyperbolic between the cusp (where Up — Vc — CgVAp/ \/ Cg + Vj) and the slow- 
magnetosonic surface (where Up — VsM,p) , elliptic between the SM and the fast-magnetosonic 
surface (where Up — Vfm,p) and hyperbolic further out. The magnetosonic phase speeds 
involved in these definitions are the usual ones, i.e. waves travelling along the poloidal field. 
Solving Eq. (1) remains a major challenge in applied mathematics: it would require to a 
priori know the locus of these surfaces, whereas they emerge as the global solution evolves. 
In practice, one either solves the time-dependent problem with full MHD codes, or uses a 
method of variable separation. 

Self-similar solutions allow to solve the full set of MHD equations without any approxi- 
mation. The problem reduces to propagate the solution along a self-similar variable x = z/r, 
which involves the inversion of a matrix. Its determinant vanishes at three singular points 
where the following numbers become equal to unity: Msm = V/VsM,n, Ma = V/VA,n and 
MpM = y/yFM,n, where V = Up ■ n and Va^u = Va • n are projections in the direction n 
(see Fig. 1) and l^jM/FM.n = 1 (KI + ± ^{Vj + Cjf - AC^Vl^) (Ferreira & Pelletier 
1995). Not all these critical points coincide with the above mentionned points where the 
flow changes type: Ma — m, Msm — Up/Vc ~ Up/VsM,p (in the cold Keplerian limit) but 
MpM < n — Up/Vfm,p- The critical FM surface is always located downstream in the hyper- 
bohc region. The necessary conditions to provide cold super-Alfvenic jets from Keplerian 
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Fig. 1. — Geometry of a disk wind configuration and definitions of quantities related to tfie 
Alfven critical surface. The unit vector n is defined as (z, r)/-\/p~+^. 

disks were given in Fcrrcira (1997). Once fulfilled, a super-A solution propagates further 
away until collimation by the hoop-stress takes place. This produces an unavoidable de- 
crease of the projected velocity {V —Ur — 0), even if the poloidal velocity reaches its 
asymptotic value of Qof'oV^X — 3. The only way to allow for a super-FM solution with 
MpM > 1 is to lower this projection effect by forcing the magnetic surfaces to remain wide 
open, namely bringing the Alfven surface closer to the disk (increase \E'a). The self-similar 
expression of the GS equation in the cold regime can be written as cos^a = R{Oa'i ^a) at 
the Alfven point (Casse & Ferreira 2000a). It is an implicit equation providing the jet 
opening angle 9 a for a given position of the Alfven surface a- In the cold limit, one has 
iasiOsM ~ tan^A(l — A~^/^) for a given magnetic lever arm A ~ and initial opening 

angle 9sm (see Fig. 1). This initial opening angle is constrained by the underlying disk ver- 
tical equilibrium. The larger angle, the larger the magnetic compression and the less mass 
is being ejected. Only angles up to ~ 45° [za ^ t a) have been proved to be possible from 
Keplerian accretion disks, either with isothermal (Ferreira 1997) or adiabatic jets (Casse & 
Ferreira 2000a), but none of these solutions can become super-FM. 

Every super-FM solutions obtained by Vlahakis et al. (2000) exhibit Alfven surfaces 
closer to the equatorial plane (i.e. ~ 60°), but being not connected to a resistive MHD 
disk, they did not have to fulfill the requirement of a quasi-static disk vertical equilibrium. 
Actually no Keplerian disk would probably survive the overwhelming magnetic compression 
imposed by the bending of the field hues (or allow the imposed mass effluvium). The 
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only possibility to have an Alfven surface closer to the disk is to break this univocal link 
between and 9sm- This implies a change of some invariants (entropy and total specific 
energy) between the disk and the Alfven surface. Physically this requires an extra force 
term in the GS equation, namely a strong outwardly directed pressure gradient in the sub-A 
region. Within the self-similar framework, it means building up a large thermal pressure, 
thus an additional heating starting above the disk. The generalized GS equation becomes 
cosOa = R{Oa] + Ri3{0a] where 

R, --^(-cosg^+ . (2) 

\gA sm(^A - Oa) 



X 



+ 




din Pa 1 
din To 7 — 1 



sin^A I Cgj^cosOA dx 



is the contribution of this additional heat flux and Pa is the ratio of the plasma pressure 
Pa = PACg ji to the poloidal magnetic pressure at the Alfven point. This equation shows 
that Pa large enough {Pa ^ 1) and Rp negative are two necessary conditions to increase 
Indeed, since 9a is always smaller than , any tendency to increase 9a leads to a lowering 
of the Alfven surface. 

At the Alfven surface, Pa — 2uj\^^ where ua — ^^ta/up^a ^ 1 (Ferreira 1997; Casse 
& Ferreira 2000a) and e = h/r is the disk aspect ratio. This general expression shows that 
any cold jet (isothermal T4 = To or adiabatic Ta -C To) always displays -C 1. In order 
to have any influence on the transverse equilibrium, this additional heating must provide 
a large increase in jet temperature, namely T4 > To/e^. The second condition {Rp < 0) 
sheds light on the required heating function. The sum of the first two terms in (2) is usually 
always negative. Indeed, energy conservation gives J^a/C^^a ^ ^prl-'- ~ t^) because of the 
tremendous cooling due to the jet expansion. Since Ta must be large, the ratio TaICI a 
is always large enough (but of order unity). The third term of (2) shows that the most 
favourable situation is the presence of additional heating mainly in the sub-A region, i.e. a 
vanishing heat flux (q • ey^ = or very small) and an already decreasing temperature (due 
to adiabatic cooling). 



3. Self-similar numerical solutions 



We follow basically the same integration procedure as in our previous works (see Casse 
& Ferreira (2000b) for more details). A heating function is assumed to be present, starting 
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at the disk surface but vanishing before the Alfven point, with an adiabatic index 7 = 5/3. 
Further up, we allow for a continuous transition to a polytropic energy equation, P oc p^. 
The X-type FM critical point allows to determine the critical value Fc of the polytropic index: 
if F < Fc thermal acceleration is too inefficient (breeze-like solution) whereas if F > Fc the 
strong decrease in enthalpy leads to a shock-like solution. Although Vlahakis et al. (2000) 
used an analogous way, our solutions strongly differ by the fact that jet invariants are fixed 
by the disk. Therefore, we need first to drastically increase the jet enthalpy before fine- 
tunning the polytropic index. As in solar wind models, we are playing around with one free 
parameter (F) whilst one should solve the full energy equation. 

Figure 2 shows a typical super-FM solution obtained with Fc = 1.45. The energy input 
required can be measured by the ratio / = J-{xa, a)/E{a) since most of the heating occurs 
in the sub-A region. Solutions displayed here required / of several 10~^, allowing to get 

— 65° with Pa — 0.1 (condition Ta ~ To/e^ is verified). Note that smaller temperature 
values would also allow super-FM jets but they would just be terminated much sooner as 
in Vlahakis et al. (2000). In general 2D flows, the "causal horizon" (here the MpM = 1 
surface) is the envelope of one of the two families of characteristics (Tsinganos et al. 1996) 
and not the surface of parabolicity n = 1 (see Fig. 3). Any perturbation occurring to the flow 
downstream to the MpM = 1 surface is unable to cross this horizon. This result is generic to 
2D solutions, the only bias introduced by self-similarity is the conical shape of such surfaces, 
not their separate existence. This has strong consequences on numerical experiments, as 
already pointed out by Ustyugova et al. (1999). To ensure the absence of feedback from 
the imposed boundary conditions, the Mach cones (defined locally as the tangents to the 
characteristics) must be directed out of the computational domain at its boundaries. 



4. Astrophysical implications 

The present computed MHD fiows are the first-ever steady-state solutions describing 
an overall accretion-ejection structure from the resistive accretion disk to the super-FM jet 
region. The strict stationarity of such accretion-ejection engines depends critically on the 
thermal properties of the sub-Alfvenic region. If the plasma pressure, measured at the Alfven 
point, is a sizeable fraction of the poloidal magnetic pressure, MHD jets from Keplerian 
accretion disks can become super-FM. In the super-FM region, the jet is always facing 
a recollimation that ends up as a shock. The further jet propagation requires numerical 
time-dependent simulations. Around a protostar, such thermal pressure gradient occurs 
whenever temperatures as high as several 10^ K are reached along the inner streamline. This 
is compatible with recent observations of blueshifted UV emission lines (Gomez de Castro 
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& Verdugo 2001) and some absorption features (Takami et al. 2002). Unfortunately, the 
heating source can only be inferred from its effects and its origin remains a crucial issue. 
For instance YSOs accretion disks are assumed to be highly magnetized so one may safely 
expect that some accretion energy is also dissipated in the upper disk layers and provides 
coronal heating (Galeev et al. 1979; Hcyvacrts & Priest 1989). This is actually shown by 
both numerical simulations (Miller & Stone 2000) and some observational indication of 
accretion powered coronae (Kwan 1997). Moreover, since the central object has a hard 
surface, the shock of the infalling material provides another source of UV radiation (as well 
as X-rays), illuminating the disk and heating the sub-A region (Ferro-Fontan & Gomez de 
Castro 2003). An alternative to this additional heating would be the presence of a high 
pressure inner flow (a "spine") forcing the MHD disk wind to open up. In YSOs such a 
flow could be provided by the interaction between the protostellar magnetosphere and the 
disk (Fcrrcira et al. 2000; Matt et al. 2002; Romanova et al. 2002). Temperatures required 
around a compact object imply a relativistic plasma. In this case, the inner pressure could 
be provided by an inner beam composed of relativistic electron-positron pairs, heated and 
accelerated inside the hollow part of the disk wind (Renaud & Henri 1998). 

On the other hand if an accretion-ejection engine cannot provide this additional heating 
or if there is no inner "spine" , then the thermal pressure is negligible at the Alfven surface 
and jets remain sub-FM. RecoUimation towards the axis leads to the formation of a shock 
and the overall structure is therefore unsteady. However, no MHD signal can propagate 
upstream along the magnetic field towards the disk. Instead, in this hyperbolic region, the 
information that a shock has occurred is first carried away by MHD waves travelling along 
the characteristics until the n = 1 surface is reached (sec Fig. 3). Once in the elliptic domain, 
the fastest mode travels along the magnetic field down to the disk. As a consequence, the 
time taken by MHD waves to inform the disk is always larger than the one given by e.g. 
computing the time r = J ds/VpM.p taken by the fast mode along the same ficldlinc. As 
an illustration, let us take the cold sub-FM solutions of Ferreira (1997), as pictured in his 
Fig. 6. For the solutions recoUimating right after the Alfven surface this time r is roughly 
equal to the orbital period To at the anchoring radius, whereas for those recoUimating much 
farther away, one gets r > IO^Tq. The presence of sporadic jet events with time scales much 
larger than disk dynamical time scales (Raga et al. (2002), Gallo et al. (2003) and references 
therein) could fit into the picture of an accretion-ejection engine trying to adjust itself. 

The amount of large scale poloidal magnetic flux trapped in accretion disks is completly 
unknown. The above astrophysical implications only hold if this flux is large enough so 
that an equipartition field spans at least one decade in radius in the disk. Indeed, in such 
circumstances, there is no physical reason for strong gradients in jets and one may expect an 
almost plane Alfven surface (Krasnopolsky et al. 1999; Casse & Keppens 2004). This kind of 
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jets display dynamical properties (acceleration and coUimation) that weakly depend on the 
radial (inner and outer) boundary conditions, as in self-similar models. On the contrary, if 
the flux is small and concentrated at the inner edge of the disk, one would expect an almost 
spherical expansion of the field lines as in X-wind models (Shu ct al. 1994). Observations 
allowing to infer jet velocity patterns and to relate them to the source (Garcia et al. 2001; 
Bacciotti et al. 2002; Pesenti et al. 2003) are necessary to discriminate between these two 
extreme pictures. 

F.C is a postdoctoral fellowship of the European Community's Human Potential Pro- 
gramme PLATON under contract HPRN-CT-2000-00153. F.C. would also like to thanks the 
team SHERPAS for its hospitahty during his stay at the LAOG. 
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Fig. 2. — Typical super-FM disk wind with ^ = 0.03, e = 0.03 {h = tr). Density, pressure 
and temperature are normalized to their value at the disk midplane, the magnetic field 
components to Bz{z = 0) and the velocities to the Keplerian speed at the anchoring radius 
Tq. All magnetic field components remain comparable from the disk surface to the Alfven 
point. Note that the density profile inside the disk, where both Uj. and are negative, is 
very different from a gaussian. RecoUimation takes place at 2; ~ 3 lO^ro. 



- 11 - 




12 3 

R (AU) 



Fig. 3. — Poloidal cross section showing the three critical surfaces (SM, A and FM), some 
characteristics (blue lines, in the hyperbolic domain) as well as the two elliptic regions 
(shaded). Contrary to Vlahakis et al. (2000), the streamlines (red lines) are computed from 
the midplane of the resistive accretion disk. This solution has ^ = 0.09, e = 0.03 and 
Fc = 1.56. 



